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ABSTRACT 


The  research  described  encompasses  five  tasks  each  of  which  treats, 
from  either  an  experimental  or  theoretical  base,  the  measurement  or 
interpretation  of  the  high- temperature  thermophysical  behavior  of  model 
*y ; terns  which  show  potential  applications  to  the  environmental  resistant 
materials  and  the  spacecraft  survivability  programs.  While  there  is 
significant  overlap  between  these  casks,  the  first  chapter  is  addressed  to 
the  environmental  resistant  materials  programs,  while  the  later  four 
chapters  are  addressed  to  the  spacecraft  survivability  program. 

Measurements  have  been  completed,  in  a  specially  designed  apparatus 
using  state  of  the  art  pvrometery,  on  the  triple  point  of  graphite  by  a 
pulse  heating  method.  The  surface  radiance  temperature  of  melting  graphite 
at  a  wavelength  of  655nm  is  4330  K,  while  the  triple  point  is  estimated  to 
be  4530  K  with  an  estimated  uncertainty  of  150  K. 

Calculational  algorithms  based  have  been  expanded  to  calculate  direct 
radiative  interactions  among  many  coupled  channels  and  provide  the 
continuum-continuum  matrix  elements  needed  in  lineshape  theory  as  well  as 
the  bound- continuum  matrix  elements  needed  to  describe  photodissociation 
phenomena.  Calculations  have  been  made  to  study  pressure  broadening  of 
atomic  lines  in  both  weak  and  strong  laser  fields.  Power  broadening  and 
saturation  of  continuum  transitions  in  intense  multiphoton  fields  have  also 
been  explored. 

A  combined  theoretical  and  experimental  program  is  being  pursued  to 
increase  the  understanding  of  the  relationship  of  structure  and  reactivity 
of  graphitic  species.  This  program  focusses  on  the  fundamental  chemical 
properties  of  very  large  polvaromatic  species.  Theoretical  methods  have 
been  developed  that  can  reliably  predict  structure  and  reactivitv  of  large 


graphitic  structures.  This  theory,  a  modified  Huckel  approach,  accurately 
predicts  enthalpies  of  formation  of  pi-radicals.  A  high  temperature,  low 
pressure  flow  reactor  has  been  built.  This  apparatus,  which  has  a 
capability  to  3000  K,  will  be  used  to  test  some  of  the  results  of  the 
theoretical  program. 

The  rates  and  mechanism  of  ionic  polymerization  in  highly 
unsaturated  organic  compounds,  structures  and  thermochemistry  of  relevant 
ions,  photodissociation  of  the  polymeric  ions,  and  the  repair  mechanisms 
which  generate  polymeric  ions  after  photodissociation  have  been  studied. 

The  work  focussed  on  the  structure  and  reactivity  of  the  phenylium  ion, 
which  is  ubiquitous  in  irradiated  aromatic  systems.  Both  the  formation  of 
this  ion  and  its  reaction  have  been  studied.  Ionic  condensation 
mechanisms,  which  occur  without  activation  barriers  have  also  been  studied. 

Development  of  apparatus  and  techniques  to  study  the  molecular  basis 
for  laser  induced  vaporization  of  refractory  materials  such  as  graphite 
continued.  The  effects  of  laser  energy,  optical  and  sample  geometry  and 
sample  morphology  on  the  distribution  of  the  vapor  phase  species  of  carbon 
ablated  from  a  selection  of  graphites  has  been  examined.  Improvements  have 
been  made  in  the  experimental  apparatus  and  the  modelling  of  the  results. 


DETERMINATION  OF  THE  TRIPLE-POINT  TEMPERATURE  OF  GRAPHITE 
FROM  MEASUREMENTS  BY  A  PULSE-HEATING  TECHNIQUE 


A.  Cezaifliyan  and  A.  P.  Miiller 


ABS^RAC  r 


Measurements  of  soectral  radiance  of  graphite  in  the  vicinity 
of  its  triple  point  have  been  performed  by  a  rapid  pulse-neating 
technique.  The  method  is  based  on:  (a)  resistively  heating  the 
specimen  in  a  pressurized  gas  environment  (-200  atm)  from  room 
temperature  to  its  melting  point  in  less  than  20  ms  by  passing  an 
electrical  current  pulse  through  it;  and  (b)  simultaneously 
measuring  the  radiance  temperature  of  the  specimen  surface  at 
655  nm  every  120  us  by  means  of  a  high-speed  pyrometer.  The 
results  for  two  different  grades  of  POCO  graphite  (DFP-1  and  AXM- 
5Q1)  are  in  good  agreement  and  yield  a  value  of  4330  +  50  K  for 
the  spectral  radiance  temperature  (at  655  nm )  at  the  triple 
point.  On  the  basis  of  the  present  value  and  literature  data  for 
normal  spectral  emittancei  the  true  triple-point  temperature  of 
graphite  is  estimated  to  be  4530  +  150  K. 


1.  INTRODUCTION 

During  the  past  decade  or  so,  there  have  been  numerous 
studies  of  the  so  1  id-1 iquid-vapor  triple-point  of  graphite 
reported  in  the  literature  (see  reviews  in  Refs. Cl, 2D).  Most 
measurements  of  triple-point  pressure  have  yielded  values  in  the 
range  of  about  100-110  atm.  However,  the  reported  measurements 
of  triple-point  temperature  are  in  considerable  d l sagreement , 
yielding  values  which  span  the  range  of  approximately  4000  to 
5000  K.  The  major  problem  in  such  measurements  arises  from  the 
high  vapor  pressure  of  graphite  at  elevated  temper atures :  as  the 
triple-point  is  approached,  a  plume  of  cooler  carbon  vapor  or 
soot  forms  near  the  specimen  surface,  thereby  severely  limiting 
the  accuracy  of  optical  pyrometry. 


A  feasibilty  study  C31  in  our  laboratory  has  shown  that 
effects  of  specimen  evaporation  on  pvrometric  temperature 
measurements  can  be  suppressed  or  at  least  minimized  under 
certain  conditions,  namely:  (1>  the  pressurized  inert  -as 
environment  of  the  specimen  must  contain  a  small  amount  of  oxygen 
to  combine  with  the  opaque  carbon  vapor,  forming  optically 
transparent  gases  CO  and  CQS;  (2)  the  specimen  heating  rate  must 
be  about  10riK-s'1  or  higher  to  minimize  the  effect  of  convective 
turbulence  on  the  sighting  path  of  the  pyrometer . 

As  a  result  of  this  study,  a  high-speed  six-wavelength 
pyrometer  was  designed  and  constructed  C4D  to  provide  a  means  of 
estimating  true  temperature  of  a  rapidly-heated  specimen  from 
accurate  measurements  of  surface  radiance.  Measurements  to 
determine  the  triple-point  temperature  of  graphite  were  then 
initiated  C5].  Unfortunately,  accurate  results  for  temperature 
were  precluded  by  a  long-term  drift  in  the  response  of  the  new 
pyrometer  discovered  during  the  course  of  successive  calibration 
experiments.  The  source  of  instability  was  subsequently  determi¬ 
ned  to  be  the  interference  filters  which  have  now  been  replaced. 

The  present  report  briefly  describes  < a )  recent  measurements 
on  two  different  grades  of  POCO  graphite  at  temperatures  in  the 
vicinity  of  the  triple  point,  (b>  results  for  the  radiance 
temperature  (at  655  nm )  of  graphite  at  its  triple  point  from 
which  an  estimate  of  the  true  triple-point  temperature  is  made. 

An  alternative  method  of  determining  true  temperature  is  also 


briefly  discussed. 


2.  MEASUREMENTS  AND  RESULTS 

Details  regarding  the  construction  and  operation  of  the 
measure, r.er  t  system,  which  includes  an  electric  power -pul  sing 
circuit?  a  h igh-pressur e  experiment  chamber,  the  multi-wavelength 
pyrometer,  and  a  digital  storage  oscilloscope,  are  given  in  the 
previous  report  C51.  A  significant  improvement  in  time  resolu¬ 
tion  (a  factor  of  six)  was  achieved  in  the  present  work  by 
replacing  the  digital  scope  with  a  new  data  acquisition  system 
capable  of  digitizing  the  pyrometer  signals  with  13-bit  resolu¬ 
tion  at  rates  of  up  to  100  kHz. 

The  first  step  in  the  two-stage  calibration  of  the  six- 
wavelength  pyrometer  has  been  completed.  It  involves  calibrating 
the  655  nm  "channel"  in  a  steady-state  experiment  with  a 
tungsten-filament  lamp  which,  in  turn,  had  been  calibrated 
against  the  NBS  Photoelectric  Pyrometer  operating  at  the  same 
effective  wavelength.  The  second  step  consists  of  transferring 
the  calibration  to  the  other  five  channels  (centered  at  500,  600, 
700,  800,  and  900  nm )  by  means  of  a  steady-state  blackbody 
furnace.  A  graphite  furnace,  which  has  been  recently  installed 
for  this  purpose,  is  currently  being  tested  to  determine  its 
blackbody  quality. 

Specimens  were  fabricated  in  the  form  of  rectangular  strips 
from  graphite  material  which  was  supplied  by  the  manufacturer  in 
the  form  of  thin  (0.5  mm)  sheets  and  designated  as  either  PQCO 
DFP-1  or  POCO  AXM-5Q1.  The  nominal  dimensions  of  each  strip  was 


25  mm  long  by  3  mm  wide.  The  center  portion  of  each  strip  was 
" °ec ked-down"  in  order  to  define  a  small  "effective"  specimen 
with  nominal  dimensions  of  2mm  long  by  1  mm  wide.  This  procedure 
-  educed  the  effective  surface  area  of  the  specimen  arq  'enc?,  the 
amount  of  carbon  vapor  released  since  melting  was  limited  to  the 
effective  portion  only. 

Prior  to  each  pulse  experiment,  a  strip  specimen  was  mounted 
vertically  inside  the  experiment  chamber  which  was  then  filled 
with  the  desired  amount  of  oxygen  (partial  pressures  in  the  range 
of  150  to  600  psig)  and  finally  pressurized  with  argon  to 
approximately  3000  psig  < ~200  atm;.  Upon  equi 1 ibrat ion  of  the 
gas  mixture,  the  circular  viewing  target  (diameter  =  0.5  mm)  of 
the  pyrometer  was  focussed  on  the  surface  of  the  "effective" 
specimen.  The  specimen  was  then  rapidly  heated  from  room 
temperature  to  its  melting  temperature  in  less  than  20  ms  by 
passing  an  electrical  current  pulse  through  it.  The  heating  rate 
varied  with  experiment  in  the  range  2x10s5  to  5x10-'  K-s-1. 

The  variation  of  surface  radiance  temperature  (at  655  nm ) 
with  elapsed  time  during  a  typical  pulse-heating  experiment  is 
presented  in  Fig.l.  The  plateau  in  the  temperature-time  function 
is  a  char ac ter i st i c  of  solid  to  liquid  transformations!  the 
glassy  appearance  of  the  remaining  portion  of  effective  specimen 
provided  further  evidence  that  melting  had  occurred.  In  all 
cases,  the  plateau  radiance  temperature,  which  is  indicated  by 
the  horizontal  line  in  Fig.l,  was  taken  as  the  maximum  temperat¬ 


ure  along  the  initial  flat  portion  of  the  plateau. 


The  results  for  plateau  radiance  temperature  determined  from 
rapid  melting  experiments  on  two  different  grades  (DFP-1  and  AXM- 
5Q 1 >  of  PQCO  graphite  are  summarized  in  Fig.  3.  As  may  ne  seen, 
the  plateau  temperature  decreases  linearly  with  oxygen  center-  in 
the  pressurizing  gas  enviroment.  The  effect  appears  to  arise 
from  the  reaction  of  carbon  in  the  specimen  surface  with  excess 
oxygen  in  the  immediate  vicinity  of  the  specimen;  the  reaction 
(combustion)  produces  additional  visible  radiation  wnich  is 
detected  by  the  pyrometer .  This  suggests  that  the  "correct" 
radiance  temperature  at  melting  may  be  determined  by  extrapolat¬ 
ing  the  measured  values  to  zero  oxygen  content. 

The  values  of  plateau  radiance  temperature  were  fitted  by  a 
1 1  near  function  of  oxygen  partial  pressure  by  means  of  the  least- 
squares  method;  the  function  is  represented  by  the  solid  straight 
line  in  Fig.  E.  Upon  extrapolation  to  zero  oxygen  content,  the 
function  yielded  a  value  of  4330  K  for  the  radiance  temperature 
(at  655  nm )  of  graphite  at  its  triple  point. 

The  major  sources  of  error  in  the  temperature  data  arise  from 
(1)  calibration  and  operation  of  the  pyrometer,  and  (E)  physical 
and  chemical  conditions  of  each  specimen. 

A  detailed  analysis  of  sources  and  magnitudes  of  the  errors 
in  temperature  measurements  with  a  similar  high-speed  pyrometer 
is  given  in  a  earlier  publication  [61.  Allowing  for  differences 
between  present  operational  conditions  and  those  in  the  earlier 
study,  we  estimate  the  maximum  error  attributable  to  pyrometry  to 


be  about  35  K . 


The  errors  arising  from  specimen  conditions  are  considerably 
more  difficult  to  assess.  A  measure  of  variability  in  specimen 
conditions  is  given  by  the  maximum  range  of  plateau  radiance 
temperatures  observed  for  a  given  oxygen  content  ' see  Fig.  2  , 

wn ich  is  about  20  K.  There  is  also  an  additional  uncertainty 
associated  with  assigning  the  present  value  of  4330  K  to  the 
triple  point  since  our  measurements  were  performed  at  -200  atm 
rather  than  at  the  triple-point  pressure  of  -100  atm.  However, 
this  error  appears  to  be  small,  less  than  2  K,  on  the  basis  of 
the  T-P  phase  for  carbon  determined  by  Bundy  [73. 

Thus,  the  total  uncertainty  (errors  due  to  pyrometry  and 
specimen  conditions)  in  the  reported  radiance  melting  temperatur 
is  estimated  to  be  not  more  thain  ^  50  K.  It  may  be  concluded 
tnat  the  radiance  temperature  (at  655  nm )  of  graphite  at  its 
triple  paint  is  4330  +  50  K. 

True  temper ature  may  be  determined  from  measured  spectral 
radiance  temperature  by  means  of  Planck’s  law  if  the  normal 
spectral  emittance  is  known  at  the  given  wavelength.  A  data 
compilation  on  thermal  radiative  properties  by  Touloukian  and 
DeWitt  C33  shows  tnat,  far  different  surface  conditions,  the 
normal  spectral  emittance  (between  650  and  665  nm )  of  graphite  l 
in  the  'ange  0.3  to  0.9  at  1500  K  and  that,  this  property  is  a 
slowly  varying  function  of  temperature.  By  e x t r apo 1  a t  i  ng  the 
trend  in  values  between  1100  and  3300  K  to  temperatures  approach 
l ng  the  triple  point,  we  estimate  the  normal  spectral  emittance 
'at  655 nm ;  of  graphite  at  its  triple  point  to  be  0.6  +  0.1.  On 


the  oasis  of  the  extrapolated  value  for  emittance  and  tne  present 
value  for  radiance  temperature*  the  true  triple-point  temperature 
of  graphite  is  estimated  to  be  4530  150  K  . 

3.  CONCLUDING  REMARKS 

Additional  pulse-heating  experiments  on  graphite  are  planned 
when  the  mu  1 t i -wave  1  eng th  pyrometer  is  fully  calibrated;  these 
must  await  final  tests  on  the  steady-state  blackbody  furnace 
which  will  be  used  to  transfer  the  calibration  of  the  a55  nm 
"channel"  to  the  other  five  wavebands*  centered  at  effective 
wavelengths  between  500  and  900  nm .  Spectral  radiance  data  in 
six  wavebands  should  yield  a  more  accurate  estimate  for  normal 
soectral  emittance  of  graphite  at  its  triple  point  and  thereby,  a 
better  estimate  for  the  true  triple-point  temperature. 

An  alternative  method  of  measuring  the  true  triple-point 
temperature  is  also  being  studied.  This  involves  similar  rapia 
melting  experiments  but  with  graphite  specimens  which  are 
fabricated  in  the  form  of  thin-walled  tubes  having  nominal 
dimensions:  25  mm  length,  3  mm  outside  diameter,  0.5  mm  wall 
thickness.  Each  tube  is  mounted  horizontally  inside  the  high- 
pressure  chamber  and  the  pyrometer  is  sighted  axial  1/  into  one  of 
the  open  ends,  thereby  approximating  b 1 ackbodv  conditions.  As 
yet,  these  experiments  have  not  yielded  meaningful  temperature 
data  due  to  problems  associated  with  carbon  vapor  condensation  at 
the  clamped  (cold  >  ends  of  che  tube.  The  results  suggest  that 
significantly  higher  heating  rates  may  be  required. 
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Vdi  ldtion  of  radiance  temperature  (at  655  nm )  of  the  specimen  sui  face 
as  a  function  of  time  during  a  typical  experiment  in  which  the 
specimen  was  pulse-heated  from  room  temperature  to  its  melting  point 
in  a  pressurized  gas  environment  (130’/.  argon*  2h'/,  oxygen)  at  JuOO  psig. 
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Broad-Band  Spectroscopy  of  Small  Molecules  in 
Intense  Laser  Fields 


Principal  Investigators:  P.S.  Julienne  and  F.  H.  Mies 


Introduction 

The  general  goal  of  this  research  is  to  provide  calculations  of 
radiative  cross-sections,  fragment  state  distributions,  and  intense  laser 
saturation  effects  for  broad-band,  wavelength  insensitive  continuum 
spectral  features  of  small  gaseous  molecules  of  possible  interest  in 
spacecraft  survivability.  Specifically  our  goal  is  to  understand  the 
interaction  of  intense  laser  radiation  with  gaseous  material  either  burned 
off  or  deliberately  released  from  a  target.  In  addition  to  providing  data 
of  possible  use  in  designing  heat  shields  and  reflective  coatings, 
knowledge  of  the  resultant  fluorescence  of  ablated  surface  material  during 
and  following  laser  illumination  of  a  target  might  serve  as  a  useful  tool 
in  discrimination  between  heavy  targets  and  lightweight  decoys. 

Theoretical  input  into  this  program  is  provided  by  extensive  sets  of 
computer  codes  available  to  the  quantum  chemistry  group.  State  of  the  art 
calculations  of  any  necessary  molecular  potentials,  couplings  and  radiative 
transition  dipoles  can  be  obtained.  Scattering  codes  provide  rigorous 
non-pertubative  radiative  scattering  matrix  elements  needed  to  evaluation 
attenuation  coefficients  for  intense,  saturating  laser  radiation 
propagating  through  gaseous  vapors.  Of  necessity,  the  program  also  involves 
extensive  analysis  of  many  new  theoretical  features  introduced  by  very 
intense  laser  fields. 

Summary  of  Accomplishments  in  FY  1986 . 

Our  original  intention  for  this  year  was  to  concentrate  on  the 
development  and  expansion  of  our  photodissociation  codes  since  we 
envisioned  this  as  the  principal  broad-band  absorption  mechanism  that  might 
contribute  to  passive  spacecraft  shielding.  However,  after  attending  the 
1985  AFOSR  conference  on  Laser-Material  Interactions  in  Albuquerque,  we 
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became  convinced  that,  because  of  the  high  pressures  and  temperatures 
associated  with  laser  ablated  vapors,  we  should  place  greater  emphasis  on 
the  study  of  the  lineshapes  of  pressured  broadened  atomic  and  molecular 
transitions  in  intense  laser  fields  as  a  possible  mechanism  for  absorbing 
incoming  laser  radiation. 

Our  scattering  programs  have  been  expanded  to  allow  for  direct 
radiative  interactions  among  many  coupled  channels,  and  now  provide  the 
tone inuum- continuum  scattering  matrix  elements  needed  in  lineshape  theorv, 
as  well  as  the  bound- continuum  matrix  elements  required  to  describe 
photodissociation .  Several  calculations  have  been  performed  to  studv  the 
pressure  broadening  of  atomic  lineshapes  in  both  weak  and  strong  laser 
fields.  In  particular  we  have  gained  extensive  experience  in  calculating 
very  accurate  nonadiabat ie  effects  in  the  conventional  weak-field  limit, 


and  have 

performed 

pre 1 im inary 

calculations 

which  exhibit 

the  power 

broadening 

and  the 

saturar ion 

of 

continuum 

transitions  in 

intense , 

multiphoton 

f ie Ids . 

Mani fescac 

ions 

of  a  .  c  .  s 

tark  shifts  and 

the  laser- 

induced  Mellow  lineshape  associated  with  the  dressing  of  fragments  by  the 
radiation  field  are  extracted  from  the  numerical  results.  This  suggests 
that  normal  perturbation  theory  as  applied  to  radiative  interactions  must 
be  carefully  scrutinized  when  applied  to  laser  attenuation  mechanisms. 

Scattering  in  the  Presence  of  Radiation 

The  general  goal  of  this  research  is  to  provide  calculations  of 
radiative  cross-sections  and  fragment  state  distributions  associated  with 
the  collisional  broadening  and  photodissociation  of  small  molecules.  This 
year  we  have  developed  extensive  experience  in  calculating  various 
nonadiabatic  effects  in  the  conventional  weak-field  limit,  and  have 
expanded  our  codes  to  treat  the  saturation  of  continuum  transitions  in 
intense,  multiphoton  fields. 


A 


2 


G3 


1 .  Non- adiabatic  Theorv  of  Atomic  Broadening  in  Weak  Fields. 

(a)  Redistribution  Calculations  for  Sr(^P«-^S)  +  Ar .  We  completed  our 
close-coupled  calculations  of  polarized  light  redistribution  in  the  Sr  +  Ar 
system.  The  calculations  predict  the  polarization  ratios  of  Sr(^P) 
fluorescence  following  line  wing  excitation  by  either  linear  or  circular 
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polarized  light.  Ab  initio  calculations  were  used  to  obtain  the  ground  and 
excited  SrAr  molecular  potentials,  which  were  adjusted  modestly  to  give 
improved  agreement  with  experiment.  The  radiative  scattering  theory  gives 
a  unified  description  of  the  absorption  coefficient  and  subsequent 
fluorescence  from  the  small  detuning  impact  limit  region  to  the  far  wings . 
The  cross  sections  for  depolarizing  collisions  of  Sr(^P)  +  Ar  were  also 
calculated.  The  calculated  absorption  coefficient,  impact  broadening  rate, 
linear  and  circular  polarization  ratios,  and  depolarization  rate 
coefficients  are  for  the  most  part  in  good  agreement  with  a  diverse  set  of 
experimental  data. 

(b)  Redistribution  Calculations  for  Na  +  Rare  Gas.  We  have  used  our 

close  -  coupled  codes  to  calculate  final  state  branching  ratios  for  fine 

2  2 

structure  states  and  alignment/orientation  for  Na(  P^S)  +  He.Ne.Ar 
systems,  and  have  obtained  good  agreement  with  experimental  measurements. 
In  analysing  the  exact  close-coupled  data  we  uncovered  a  useful  spin¬ 
uncoupling  approximation  of  the  collision  dynamics  and  its  dependence  on 
collision  velocity. 

(c)  Lorentzian  Pressure  -  Broadened  widths  for  Na.  We  have  calculated 

the  Lorentzian  pressure-broadened  widths  for  both  the  D1  and  D2  lines  of  Na 

broadened  by  collisions  with  He,  Ne ,  and  Ar .  In  addition  we  have  performed 

the  first  calculation  of  the  so-called  asymmetry  parameter  which  has  been 

widely  observed  in  the  impact  regions  of  atomic  lineshapes.  The  absorption 

2 

profile  within  the  impact  limit  for  one  of  the  isolated  P.  fine  structure 
lines  may  be  written  as 


I  ( A  ) 


7 j  (hj  )A 


(A.-5.)2  +  y2(0) 
J  2  2 


where  A^  and  5  are  the  respective  detunings  and  shift  of  the  line.  The 
detuning -dependent  quantity  7  (A^ )  may  be  shown  to  have  a  finite  limit  as 

A.-O  and  exactly  equals  the  usual  impact  limit  broadening  coefficient, 

j 

including  the  effects  of  both  elastic  and  inelastic  collisions.  The 
dependence  on  detuning  gives  rise  to  the  Lorentzian  dispersion  asymmetry 
which  has  been  observed  for  Na  +  rare  gas  systems.  We  have  calculated  the 
radiative  scattering  cross  section  for  detunings  near  the  D1  and  D2  lines 
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for  Na  broadened  by  He,  Ne  and  Ar .  The  widths  and  asymmetry  parameters  are 
qualitatively  in  agreement  with  measured  values.  Deficiencies  in  the 
existing  molecular  potentials  limit  the  expected  accuracy  of  the  scattering 
calculations . 

2 .  Nonadiabatic  Theory  of  Redistribution  in  Strong  Fields 

(a)  Collisional  Redistribution  of  Resonant  Radiation.  Close-coupling 
codes  developed  for  scattering  in  the  presence  of  radiation  fields  have 
been  used  to  calculate  the  pressure  broadening  and  collisional 
depolarization  of  the  Sr  resonance  transition  in  strong  saturating  fields 

Sr(1S)  +  Ar  +  -  Sr*^)  +  Ar  -  Sr(1S)  +  Ar  +  fiu 

As  we  increase  the  Rabi  frequency  and  decrease  the  detuning  A  of  the 
incident  laser  three  observable  properties  associated  with  optical 
collision  processes  are  obtained  from  such  calculations: 

(i)  Impact  Limit  Lineshape .  As  A^-»  0  we  can  demonstrate  that  proper 
pressure  broadened  half  widths,  which  include  both  elastic  and  inelastic 
collisional  effects,  can  be  rigorously  extracted  from  close-coupled 

I 

radiative  scattering  matrix  element.  The  first  manifestations  of  a.c.  stark 
shifts  associated  with  the  dressing  of  fragments  by  the  radiation  field  are 
extracted  from  the  numerical  results, and  confirms  the  anticipated  power 
broadening  of  the  Lorenztian  absorption  profile.  This  suggests  that  normal 
perturbation  theory  as  applied  to  radiative  interactions  must  be  carefully 
scrutinized . 

(ii)  Inelastic  transitions  involving  dressed  atomic  states.  As  increases 
and  eventually  exceeds  A^  the  asymptotic  Sr  atomic  states  become  dressed. 
Our  fully  quantal  inelastic  cross-sections  can  be  used  to  calculate 
population  dedistribution  rates  in  the  presence  of  intense  fields. 

(iii)  Collisional  redistribution  spectra  in  saturating  fields.  In  strong 
fields  we  obtain  the  redistribution  spectra  for  fiu>  in  the  impact  limit  and 
find  the  typical  three-peaked  Hollow  lineshape  due  to  fluorescence  between 
dressed  atomic  states.  Both  the  half -width  and  the  polarization  of  the 
emission  peaks  are  calculated. 

(b)  Analysis  of  Hultiphoton  Lineshape  Theory  in  Intense  Fields.  One 
of  the  most  critical  ingredients  in  understanding  the  interaction  of 
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colliding  atoms  and  molecules  with  intense,  multimoded  laser  fields  is  the 
calculation  of  collisional  lineshapes  and  the  associated  relaxation  rates 
for  the  optically  pumped  energy  levels.  In  particular  these  are  needed  to 
define  the  redistribution  spectra  which  might  serve  as  a  diagnostic  of  the 
laser- induced  vapors.  Two  complimentary  theoretical  approaches  are  used 
to  describe  these  phenomena.  The  first  is  represented  by  our  "optical 
collision"  approach  here  at  NBS  where  we  calculate  the  radiative 
scattering  matrix  elements  for  atom-atom  scattering  in  the  presence  of 
intense  monochromatic  laser  fields.  The  calculations  incorporate  manv 
subtle  and  sophisticated  non- adiabatic  and  inelastic  effects  which  are 
often  neglected  in  conventional  lineshape  theory.  However,  these  data  that 
we  generate  must  then  be  incorporated  into  very  subtle  and  complicated 
theoretical  expressions  for  the  density  matrix  of  the  interacting  molecules 
in  the  field  of  the  pumping  laser  fields.  Only  then  do  our  calculations 
yield  the  proper  description  of  the  redistributed  radiation. 

The  second  theoretical  approach  is  exemplied  by  the  tetradic 
formalism  developed  by  Prof.  Ben-Reuven  at  Tel  Aviv  University  using  a 
Liouvillian  analysis  of  the  molecular  density  matrix  in  the  presence  of 
collisions.  Fortunately,  the  formal  ingredients  of  this  sophisticated 
theory  are  the  very  scattering  matrix  elements  that  are  generated  by  our 
optical  collision  approach. 

We  have  begun  a  collaboration  with  Prof.  Ben-Reuven  to  explore  the 
relationship  between  the  tetradic  description  of  redistribution  and  the 
close-coupled  theory  of  scattering  in  strong  multiphoton  fields.  So  far  we 
have  successed  in  showing  the  explicit  relationships  required  to  analyze 
resonance  fluorescence  and  resonant  Raman  lineshapes. 

Atomic  and  Molecular  Scattering  Theory 

1.  Half -Collision  Amplitudes  -  Analysis  of  Exact  Close - Cour led  Codes.  We 
have  previously  used  a  multichannel  quantum  defect  analysis  ^MCQDA)  of 
exact  close-coupled  radiative  matrix  elements  to  show  that  these  can  often 


be  factored  as  follows 


I.-- 


where  the  'incoming'  half -collision  matrix  N...  tracks  the  non- adiabatic 

11 

scattering  of  the  incident  collision  complex,  and  summarizes  the 

inelastic  behavior  of  the  'outgoing',  radiatively  excited  collision 

complex.  We  have  speculated  on  the  unitaritv  of  these  matrices  which  would 

justify  the  persistent  used  of  the  pure  adiabatic  Born-Oppenheimer  matrix 
BO 

element  S^,^,  in  evaluating  total  radiative  cross-sections  for  continuum 
transitions.  Based  on  a  convincing  but  admittedly  indirect  numerical 
analysis  of  several  previous  studies  we  have  been  very  confident  in 
advancing  the  merits  of  this  approximation.  We  have  now  modified  our  exact 
close-coupled  codes  to  extract  the  exact  half  collision  amplitudes  without 
any  commitment  to  unitar ity.  This  new  general  capability  of  obtaining  N+ 
and  N  directly  from  close -coupled  codes  will  enjoy  a  wide  range  of 
applicability  in  continuum  spectroscopy  and  allow  us  to  judge  the  validity 
of  the  unitaritv  condition. 
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2.  Multichannel  Semiclassical  Half  Collision  Amplitudes.  In  addition  to 
solving  the  close  -  coupled  equations  exactly,  we  have  developed  a  new 
simplified  numerical  procedure  designed  to  obtain  the  half  collision 
amplitude  whenever  the  nov-adiabatic  couplings  are  restricted  to 
classically  accessible  regions  and  the  unitarity  of  N+  is  insured.  This  is 
equivalent  to  solving  a  set  of  coupled  semiclassical  equations.  We  begin 
by  choosing  a  set  of  diagonal  reference  potentials  tailored  to  the  nature 
of  the  exact  interaction  matrix.  These  are  used  to  generate  a  set  of 
elastically  scattered  WKB  reference  wave funct ions .  Two  obvious  choices 
that  we  have  explored  in  some  detail  are  either  pure  diabatic ,  or  pure 
adiabatic  reference  potentials,  although  in  any  given  situation  the 
'preferred'  choice  is  usually  self  evident.  By  using  a  random  phase 
approximation  to  neglect  quickly  oscillating  contributions  to  the  close 
coupled  equations  we  can  develop  a  set  of  first  order  differential 
equations  for  the  half  collision  amplitudes  which  are  perfectly  unitarv, 
and.  in  appropriate  and  well  understood  circumstances,  reproduce  the  exact 
close-coupled  amplitudes.  This  multichannel  semiclassical  approximation  is 
presently  being  applied  to  a  number  of  interest  situations,  including 
multichannel  curve  crossings,  and  asymptotic  frame  transformation  among 
fine-structure  states. 
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3.  Multichannel  Photodissociation  amplitudes.  The  expansion  of  our 
multichannel  close-coupling  codes  to  calculate  photodissociation  cross- 
sections  is  nearly  complete.  This  gives  very  stable  procedures  for 
evaluating  bound- continuum  transitions,  including  all  possible  subtleties 
associated  with  non-adiabatic  effects  in  molecular  dissociation.  A 
renormalized  Numerov  algorithm  has  been  implemented  into  the  codes  for  the 
calculation  of  the  usual  continuum-continuum  S-matrix  elements  needed  for 
coiiisional  problems,  but  more  specifically  for  calculation  of  bound- 
continuum  transition  matrix  elements  needed  for  photodissociation.  This 
algorithm  is  still  being  tested.  In  particular  the  code  is  designed  to 
calculate  the  final  state  distribution  of  the  photofragments.  Ultimately 
these  exact  matrix  elements  can  be  compared  to  the  predictions  of  our  half 
collision  amplitudes  based  on  MCQDA. 


Papers  Completed  in  FY86 


Three  papers  partially  supported  by  this  research  contract  have  been 
accepted  or  submitted  for  publication. 

(1)  "Nonadiabatic  Theory  of  Atomic  Line  Broadening:  Redistribution 

'  11 

Calculations  for  Sr(  P«-  S)  +Ar."  P.S.  Julienne  and  F.  H.  Mies,  to  be 
published  in  Physical  Review  A. 

(2)  "Radiative  Scattering  Close  Coupled  Calculation  of  Lorentzian  Width 
and  Asymmetry  Parameters  for  Na  +  Rare  Gas  Systems"  L.  Vahala  and  P.S. 
Julienne,  submitted  to  Physical  Review  A 

(3)  "Application  of  Recoil  and  Frame  Transformation  Models  for  Coiiisional 
Redistribution  of  Light"  P.S.  Julienne,  L.  Vahala  and  M.D.  Havey, 
submitted  to  Physical  Review  A 


GROWTH  AND  DECOMPOSITION  OF  GRAPHITIC  STRUCTURES 


Principal  Investigators:  S.  E.  Stein  and  R.  L.  Brown 

Introduction 

Graphitic  materials  are  the  most  commonly  discussed  laser-armor  protective 
materials  for  spacecraft.  Upon  laser  heating,  the  molecular  structures  in  these 
materials  decompose  and  reform,  probably  in  both  the  vapor  and  condensed  phases. 
At  present  there  exists  no  molecular  understanding  of  these  chemical  processes. 
This  program  intends  to  develop  this  understanding  in  order  to  assist  both  the 
rational  analysis  of  laser- induced  ablation  processes  as  well  as  the  development 
of  improved  protective  materials. 

The  nature  and  reactivity  of  the  chemical  structures  formed  in  these 
processes  can  control  the  effectiveness  of  laser  protective  materials  in  several 
ways.  (1)  The  formation  rate  and  nature  of  the  volatiles  produced  may  depend  on 
both  the  nature  and  concentrations  of  precursor  sites  in  the  solid  as  well  as  on 
the  reaction  mechanism.  (2)  The  chemical  structure  and  reactivity  of  the 
materials  remaining  after  ablation  by  the  initial  laser  pulse  will  be  determined 
largely  by  chemical  factors.  (3)  The  condensation  of  carbon  vapors  to  form 
refractory  particles  may  not  be  controlled  by  diffusion,  but  by  chemical 
reaction,  since  reformed  structures  must  form  strong  chemical  links  to  resist  re 
evaporation . 

The  development  of  realistic  working  models  for  the  above  reactions 
requires  a  substantial  increase  in  our  understanding  of  the  structure  and 
reactivity  of  graphitic  species.  To  develop  this  understanding  we  propose  to 
continue  our  combined  theoretical/experimental  program  which  focusses  on  the 
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fundamental  chemical  properties  of  very  large  polyaromatic  structures.  Graphite 
and  nearly  all  carbonized  organic  materials  are  composed  primarily  of  these 
structures . 


Early  in  this  program  we  discovered  that  the  atomic  orientation  at  the  edgei 
of  these  structures  and  their  overall  size  were  critically  important  factors  in 
determining  the  reactivity  of  these  structures.  Our  theoretical  work  has 
therefore  been  focussed  on  separating  and  examine  individually  these  two  vari¬ 
ables.  Our  experimental  program  is  devoted  to  examining  well  defined  carbon 
deDosition  reactions  which  can  be  connected  with  our  theoretical  work. 
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Results 


Theory 

Our  major  finding  during  the  last  year  has  been  that  a  modified  version  of 
classical  Huckel  theory  is  both  qualitatively  reliable  and  quantitatively 
accurate  when  applied  to  large  polyaromatic  species.  In  addition,  we  found  that 
both  the  simpler  theories  (PMO,  structure-resonance  theory  and  other  perturbation 
methods)  and  the  more  computationally  involved  SCF-methods  make  serious  errors  in 
their  extension  to  very  large  molecules.  We  now  have  a  predictive  theory  that  we 
feel  is  not  prone  to  major  errors.  Fortunately,  this  theory  is  sufficiently 
simple  to  allow  calculations  for  molecules  as  large  as  we  wish  (up  to  several 
thousand  carbon  atoms) .  One  important  result  of  this  application  is  that  attack 
of  interior  carbon  atoms  in  graphitic  sheets  results  in  a  significant  loss  of 
resonance  energy  and  is  thus  not  likely  to  be  an  important  process  in  graphitic 
systems.  Another  finding  is  that  large  differences  in  reactivity  between 
different  edge  sites  remains  even  in  the  limit  of  an  arbitrarily  large  graphite 
structure.  Simpler  theories  incorrectly  predicted  that  all  sites  become  equally 
reactive  in  the  limit  of  very  large  cluster  size.  More  sophisticated  theories 
incorrectly  predict  a  finite,  edge-dependent  band  gap  in  very  large  clusters  as 
well  as  an  unreasonable  bond  alteration.  We  have  also  discovered  that  this 
modified  Huckel  method  accurately  predicts  heats  of  formation  of  pi-radicals. 

This  technique  has  not  previously  been  thoroughly  applied  to  these  species,  so 
that  this  observed  accuracy  strongly  supports  the  reliability  of  this  theoretical 


approach . 


In  general  terms,  edge  carbon  atoms  behave  very  much  like  those  of  normal 
polynuclear  aromatic  molecules,  with  properties  depending  strongly  on  the  loca_ 
structure.  Phenanthrene- like  edges  are  rather  unreactive  and  thermodynamically 
stable.  They  are  as  susceptible  to  free  radical  and  electrophilic  attack  as  are 
structurally  related  positions  in  naphthalene  and  phenanthrene .  Anthracene  - 1 ik. 
edges  are  thermodynanamically  unstable  and  highly  reactive,  with  properties  coir- 
parable  to  those  of  the  linear  polyacenes.  Regardless  of  edge  structure,  in¬ 
terior  carbon  atoms  situated  only  several  bond  lengths  from  an  edge  have  pro¬ 
perties  (bond  orders,  free  valence)  not  far  from  those  of  the  interior  carbon 
atoms  in  cata-condensed  polyaromatic  molecules  like  coronene .  A  significant 
general  finding  for  the  molecules  studied  is  that  there  is  only  a  weak  relation 
between  the  position  of  the  highest  occupied  molecular  orbital  and  the  reactivity 
of  the  most  reactive  position.  For  instance,  while  very  large  Series  3 
molecules  (phenanthrene- like  edges,  see  Fig  1)  will  have  higher  HOMO  energies 
than  small  Series  1  molecules  (anthracene- like  edges,  see  Fig  1)  ,  they  will  be 
far  less  reactive  (cf.  ,  Figs.  2  and  3)  .  This  HMO  method  predicts  that  thermody¬ 
namic  stability  is  strongly  influenced  by  edge  type.  In  contrast,  because  they 
utilize  local  properties  only,  group  additivity  methods  do  not  predict  signifi¬ 
cant  edge  effects. 

Wave  functions,  which  describe  overall  molecular  properties,  separate  into 
two  distinct  classes  for  molecules  with  reactive  edges.  Those  of  a  more  stable 
variety  are  qualitatively  similar  to  those  of  an  infinite  plane  and  of  the 
molecules  with  stable  edges.  Those  of  a  less  stable  variety  correspond  to 
orbitals  whose  electron  density  are  located  primarily  at  the  edges.  The  latter 
variety  lead  to  a  near  zero  band  gap,  and  the  associated  spectroscopic  and 
electronic  properties,  for  molecules  of  relatively  small  size.  For  example,  a 


molecule  with  reactive  edges  containing  60  carbon  atoms  has  the  same  band- gap  as 
a  molecule  of  the  stable  variety  containing  between  2000  and  2500  carbon  atoms. 

Experiment 

Most  of  the  year  was  devoted  to  the  construction  of  a  very  high  temperature, 
low  pressure  flow  reactor.  A  commercial  graphite  element  furnace,  capable  of 
achieving  temperatures  up  to  3000  K  has  been  modified  and  connected  to  a 
specially-designed,  differentially  pumped  mass  spectrometric  detection  system. 
Either  molecular  beam  sampling  or  indirect  detection  is  possible.  The  reactor  is 
made  of  high  density  graphite,  with  a  volume  of  ca.  200  cm^  and  an  exit  aperture 
of  3  mm  diameter.  This  relatively  large  size  will  enable  operation  over  a  wide 
range  of  well  controlled  pressure  and  chemical  conditions.  The  inlet  and  flow 
calibration  system  is  near  completion  and  will  allow  independent  flows  of  two 
reactants.  Temperature  measurement  can  be  made  by  either  thermocouples  or 


optical  pyrometry. 


Figure  Captions 


Figure  1.  Homologous  series  of  hexagonally  symmetric  molecules. 


Figure  2.  Maximum  free  valence  for  each  molecule  in  the  four  homologous  series 
as  a  function  of  molecular  size. 


Figure  3.  Energy  of  highest  occupied  molecular  orbital  (HOMO)  as  a  function  of 
molecular  size.  The  reciprocal  of  the  square  root  of  the  number  N  of 
pi-electrons  in  the  molecule  was  used  for  the  abscissa.  This  was  done  to  expand 
the  scale  near  the  origin  (i.e.  ,  where  N  goes  to  infinity)  .  The  value  of  the 
Huckel  electron  repulsion  integral  is  the  energy  unit. 
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Figure  2.  Maximum  free  valence  for  each  molecule  in  the  four 
homologous  series  as  a  function  of  molecular  size. 
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APPLICATION  TO  SPACECRAFT  SURVIVABILITY 


Sharon  G.  Lias 


For  understanding  the  generation  and  dispersal  of  smokes  and 
jjstur  in  4,  materials  in  space,  the  kinetics  of  evaporation  of  refrac¬ 
tory  materials  under  laser  irradiation,  and  condensation  processes 
leading  to  the  formation  of  microscopic  particles,  graphitic  materials 
and  polyarylacetylenic  polymers  are  of  particular  interest.  This 
project  involves  the  molecular-level  study  of  ionic  polymerization 
median  isms  in  polyunsaturated  linear  and  aromatic  compounds,  the 
structures  and  thermochemistry  of  relevant  ions,  the  photodissociation 
or  tae  polymeric  ions,  and  the  repair  mechanisms  which  regenerated  the 
polymeric  ions  after  photodissociation.  The  model  systems  studied 
previously  (diacetylene,  cyanoacetylene,  cyanogen,  phenyl  acetylene  and 
styrene)  ’  ’  were  focussed  on  highly  unsaturated  organic  compounds  in 
which  positive  ions  are  generated  and  undergo  rapid  polymerization 
with  surrounding  neutral  molecules  (which  can  be  identified  with 
molecules  vaporized  from  a  graphitic  surface  material).  The  product 
ions  from  these  ion/ molecule  reactions  are  larger  unsaturated  species 
which  nave  broader  absorption  spectra  and  thus  are  more  desirable  as 
spacecraft  obscurants. 


Conclusions  derived  from  earlier  work  of  this  project  are  that 

+ 


the  highly  unsaturated  fragment  ions  (for  example,  C  ^  Hi  and  C^H->  in 
phenylacetylene)  react  at  every  collision  with  the  precursor  aromatic 
molecules,  and  tnat  efficient  condensation  reactions  occur  followed  by 
elimination  of  a  molecule  of  acetyene.  The  parent  radical  cations  in 
these  systems  react  with  the  corresponding  neutral  molecules  to  form 
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condensation  L  0  n  s  ,  which  then  dissociate  by  the  elimination  of  a 
stable  molecule  (benzene,  hydrogen,  ethylene,  or  acetylene  for  ex¬ 
ample). 

The  ICR  cell  is  fitted  with  a  quartz  window,  and  the  ions  can  be 
irradiated  with  (a)  a  xenon  arc  lamp  fitted  with  various  interference 
filters  which  allow  one  to  irradiate  the  trapped  ions  with  monochroma¬ 
tic  light,  and  (b)  an  argon  ion  laser.  Studies  on  the  photodissocia¬ 
tion  mechanisms  of  parent  radical  cations  and  ionic  products  of 
ion, molecule  reactions  have  been  carried  out.  In  these  studies  three 
different  types  of  behavior  have  been  seen:  1)  the  loss  of  small 
molecules  which  shorten  the  carbon  chain,  2)  the  loss  of  a  molecule  of 
hydrogen  resulting  in  an  increase  in  the  degree  of  unsaturation,  3) 
the  loss  of  monomer  molecules  which  is  equivalent  to  depolymerization. 
Further,  the  photodissociation  mechanisms  of  the  ions  have  been  direc¬ 
tly  determined,  and  the  modes  of  reaction  and  structures  of  dissocia¬ 
tion  products  have  been  examined.  At  the  present  time,  little  is 
known  about  photodissociation  mechanisms  of  ions  formed  as  products  of 
ion/' molecule  reactions. 

Work  during  FY8b  focussed  on  two  problem  areas:  (a)  the  chemis¬ 
try  of  tne  phenyliura  ion  (c-C^H^*),  which  is  ubiquitous  in  irradiated 
aromatic  systems;  and  (b)  on  detailed  examinations  of  the  actual  mech¬ 
anisms  of  ionic  condensation  processes. 

Structure  and  Reactivity  of  the  Phenyliura  Ion 

Until  our  recently  completed  work.4,  little  was  known  ibout  the 
structure  and  reactivity  of  the  phenyliura  ion.  Our  results  have  demon¬ 
strated  that  previous  conclusions  derived  from  c oil  i  si  on-in d need - 
dissociation  ( C I L) )  experiments  were  incorrect.  That  is,  earlier  re- 


suits  hud  led  to  the  conclusion  that  acyclic  C  b  H  ^  +  ions  are  reactive 
with  acetylene,  and  phenylium  ions  are  unreactive;  our  study  proves 
conclusively  that  the  cyclic  are  reactive,  while  the  acyclic 

isomers  are  unreactive  with  aromatic  molecules  and  aikynes.  Rate 
constants  for  reactions  of  CbH^  +  had  not  previously  been  measured. 

From  ab  initio  calculations,  at  was  known  that  the  most  stable 
1 o  d  • +  ion  is  the  singlet  phenylium  i^n,  while  acyclic  isomers  are 
predicted  to  have  heats  of  formation  at  Least  15  kcal/moi  higher. 
From  tnis  knowledge  of  the  energetics,  it  necessarily  follows  that 
phenylium  ions  are  uniquely  produced  in  the  dissociative  protonation 
of  n  a  1  o  benzenes: 

MH+  +  C  b  H  5  X  ->  c-Cbrt5+  +  IIX  +  M  (1) 

Fnenyiium  ions  must  also  be  uniquely  produced  in  the  threshold  elect¬ 
ron-impact  dissociation  of  the  halo benzenes: 

CbH5X  +  e  ->  c-C6d5+  +  X  +  2e  (2) 

or  in  tue  Low  energy  photodissociation  of  the  halo  benzene  cations: 

C6H5X+  +  hv  _>  c"C6H5+  +  X  (3) 

That  at  Least  two  isomeric  Cb4^  +  ions  are  formed  by  electron 

impact  fragmentation  of  organic  compounds  as  well  as  in  ion /molecule 
reactions  is  demonstrated  by  the  results  given  in  Figure  1.  Figure  1, 
which  represents  the  time -dependent  decay  of  ions  formed  in  the 

e . e c t r o n - i m p a c t  fragmentation  of  C^H^Cl,  clearly  shows  that  CbHb  + 
po  p  u  L  it  ions  or  distinctly  different  reactivities  exist.  The  reactive 
2^4  iins  exhibit  a  rate  constant  of  5  x  l')~  1  J  for  the  reaction: 

o  b  4  b  +  +•  C  b  4  ^  C  1  ->  C  i  2  H  i  ( ^  +  Cl  (4) 

The  C  b  4  b  ions  produced  in  the  reaction: 

C  ^  H  j  2  ^  4  3  4  r  -  >  C  h  4  ^  H  B  r  (  5  ) 
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exhioit  the  same  race  constant  for  reaction  with  chlorobenzene,  and 
therefore  can  be  assumed  to  have  the  same  structure  as  the  reactive 
population  of  fragment  ions  in  the  aromatic  system. 

In  order  to  assess  which  of  the  two  populations  corresponds  to 

the  phenylium  ion,  non-fragmenting  C^H^C1+  ions  produced  with  11  eV 

electrons  were  irradiated  witn  2  eV  photons  of  tne  argon  ion  laser  for 

a  short  period  of  time,  after  which  their  decay  rate  was  monitored. 

Under  tnese  conditions,  only  one  population  of  ions  is  produced; 

criose  o  H,"*’  ions  react  with  the  parent  chlorobenzene  molecules  with  a 
oo  r 

reaction  rate  constant  (5  x  10  ^  cm^/molecule-s)  which  is  identical 
to  that  obtained  for  the  more  reactive  population  formed  in  the  elect¬ 
ron  impact  fragmentation  experiment  (Figure  1).  The  photodissociation, 
which  occurs  through  a  two-photon  process,  can  only  lead  to  the  forma¬ 
tion  of  phenylium  ions,  and  therefore,  it  can  be  concluded  that  the 
reactive  population  has  the  c  -  C  ^  H  3  +  structure. 

Because  of  the  substantial  differences  in  energy  and  energy 
barriers  between  the  phenylium  ion  and  the  acyclic  isomers,  the  iso¬ 
meric  distribution  of  species  formed  under  any  given  conditions 

will  tie  determined  by  the  energy  content  with  which  the  ions  are 
formed.  Thus,  in  acetylene,  where  C  ^  H  3  +  ions  are  generated  in  a  simple 
addition  reaction: 

C  3  +  +  C2H2  ->  C6H5+  (6) 

c.ie  product  will  originate  with  a  total  energy  of  3  5  3  kcal/mol; 

cats  is  d2  xcai/mol  above  the  ground  state  of  the  phenylium  ion,  but 
is  not  sufficient  to  cause  dissociation.  Reaction  5,  which  is  the 
main  secondary  reaction  channel  in  acetylene,  was  found  to  yield  two 
C  d-4"  populations  (a  reactive  and  an  unreactive),  with  the  unreactive 
c  H  +  species  accounting  for  60%  of  the  total  population.  although 


30 


CIO  experiments  had  led  to  the  conclusion  that  the  latter  ions  are  the 
pnenyliura  ions,  authentic  phenyl  ions  generated  through  reaction  4 
reacted  with  acetylene  in  a  fast  reaction.  The  rate  constant  for  this 
reaction  in  pure  acetylene  was  identical  to  that  determined  tor  the 
authentic  phenyl  ions,  thus  demonstrating  conclusively  that  the  react¬ 
ive,  rather  than  the  unreactive,  population  has  the  phenylium  struc¬ 
ture.  The  rate  constant  which  can  be  ascribed  to  the  reaction  of  the 
aevclic  ions  with  acetylene  is  0.15+0.02  x  10  ^  cra^/molecule-s, 

-  0  J  — 

wnich  is  within  experimental  error  the  same  as  that  obtained  for  the 
unreactive  Cbri3  +  ions  (Figure  1)  in  C^H^Cl  when  C,D,  is  added  to  the 
system. 

In  the  reaction: 

C  3  H  j +  +  C3H3C1  ->  C6H5+  +  HC1  (7) 

the  product  ions  are  formed  with  somewhat  less  excess  energy  over  that 
if  the  ground  state  phenylium  ion  (<7  2  kcal/'iol)  ;  any  excess  vibra¬ 
tional  energy  imparted  to  HC1  will  reduce  the  energy  content. 

The  relatively  high  abundance  of  the  unreactive  population  indicates 
that  most  of  the  energy  is  retained  by  the  C  ^  H  3  product  ton. 

Ionic  Condensation  Mechanisms 

Another  study"*  explored  the  mechanisms  of  ionic  condensation 
processes  which  occur  without  activation  barriers  in  competition  with 
t  is:  proton  transfer  reactions: 

/  —  >  3  H +  +  A  (da) 

+  / 

AH  +  S 

—  >  1  A  B  H  J  +  (3b) 

The  conclusions,  in  agreement  with  the  results  of  an  earlier  study 
rrom  this  laboratory^5,  were  that  the  relative  probabilities  or  the  two 


channels  were  determined  by  the  initial  configurations  of  the  react¬ 
ants  in  the  collision  complex,  i.e.  the  reactions  involve  two 
independent  transition  states.  Evidence  is  also  presented  that  when 
the  reactant  ion  has  excess  kinetic  energy,  a  long-lived  condensation 
ion  is  not  formed,  but  rather  a  mechanism  involving  a  "direct"  diss- 
ociiclon  occurs* 
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Time  dependence  of  Che  abundance  of  C ^ H 5  formed  in 
chlorobenzene  at  a  total  pressure  of  1  x  10”®  torr  , 
electron  energy,  40  eV.  From  a  comparison  of  the 
abundance  of  the  unreactive  population  seen  at  long  time 
with  the  total  initial  abundance  of  ions,  it  can  be 
estimated  that  approximately  14 Z  of  the  ions  are  in  the 
unreactive  configuration. 
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MOLECULAR  BASIS  FOR  LASER- INDUCED  VAPORIZATION  OF  REFRACTORY  MATERIALS 


J.W.  Hastie,  D.W.  BonneLI,  and  ?.K.  Schenck 


Ceramics  Division 


National  Bureau  of  Standards 


Introduction 


The  primarv  goal  of  this  task  is  the  development  of  a  mo lecular - ieve 1 
understanding  of  the  vaporization  and  related  ionization  and  aggregation 
processes  which  occur  when  high-power  laser  energy  interacts  with 
refractories.  Materials  of  current  interest  include  graphite,  and 
refractor;.’  carbide  materials  which  are  candidate  protective  material's  for 
orbital  applications.  The  task  rationale  and  approach  has  been  discussed  in 
detail  in  an  earlier  report,  NBSIR  8*4-2983  .  and  recent  results  are 
described  in  the  paper  "Laser  Vaporization  Mass  Spectrometry",  submitted  to 
Nucl.  Instr.  and  Methods  in  Phvs .  Res.  B.  Beam  Interactions  with  Materials 
and  Atoms . 


The  basic  experimental  approach  consists  of  coupling  a  high-power 
Nd:YAG  laser  with  a  specially  designed  very  high  pressure  sampling 
molecular  beam  mass  spectrometer  system  17HPMS ) .  The  vacuum  svstem  provides 
for  extremely  high  pumping  speed,  simulating  near-earth  orbit  gas  dvnami.cs 
even  when  laser-evolved  get.  plumes  from  refractor/  surfaces  reach  pressures 
above  one  atmosphere.  Using  laser - svnchronized  detection  methods,  time  and 


mass - reso Ived  spectra  of  lase r - induced  vapor  species  are  obtained.  The  mass 
spec  trome  tr  ic  ''MS)  data  provide  the  species  specific  i  tv  needed  for 
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a  r. 


determination  of  kinetic  and  thermodynamic  parameters  of  the  Laser  heating 


i  9 


and  vaporization  process  when  coupLed  with  temperature  and  gas  dynami' 


.nr  or mat  ton . 


if  tar"  icul  -.r  current  'interest  are  the  etfec*“s  or  laser  -tar;" 
optical  and  sampling  geometry,  and  surface  morphology  on  the  vapor  species 
distributions  .with  respect  to  concentration  and  time)  for  a  selection  of 
graphites.  Since  surface  temperatures  in  excess  of  4000  K,  and 
heating,  cooling  times  of  milliseconds,  or  less,  are  involved,  new,  or 
improved,  measurement  methods  are  required.  Methods  under  development 
include  techniques  for  determining  surface  and  plume  temperatures  and  the 
time -his tor;.’  of  vapor  species  evolution. 


Results  for  FY-86 


Modification  of  Plume  Sampling  Geometry 


Our  previous  experimental  work  confirmed  equilibrium  C  distributions 
in  the  vapor  plume  for  off-axis  (>  15°  from  surface  normal)  beam-sampling 


s 
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geometries.  Also,  surface  temperature  appeared  to  be  independent  of  Laser 
energy  at  surface  energy  Loadings  of  10  -  40  mJ/cm2 .  The  primary  observed 


effect  of  increased  laser  power  was  a  Linear  increase  in  species 
intensities  but  with  no  indication  of  temperature  changes  in  the  time-of- 
flight  or  species  distribution  data.  This  result  suggests  that  a 
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phys iochemical  restraint  is  limiting  the  temperature  increase  of  the 
surface.  Qualitative  observations  of  deposit  patterns  indicated  that  the 


beam  was  highly  forward  collimated.  Since  the  existing  apparatus  geometry 
prevented  pe  rper.dicuiar  i  beam  -  samp  i  ing  axis  to  pLume  axis  angle  of  zero 
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degrees)  sampling,  and  there  was  evidence  of  laser-plume  interaction  even 
at  relatively  low  laser  attack  angles,  the  interaction  was  modified  to 
allow  sampling  at  angles  with  respect  to  the  surface  normal  of  between  ~0° 
and  153 .  The  chamber  was  re-constructed  with  multiple  access  ports  and  an 
internal  beam  periscope  but  the  off-axis  laser  impingement  geometry  was 
retained  to  minimize  laser-vapor  interactions  and  permit  sample  rastering 
with  the  in-vacuum  two-axis  translation  stages.  The  new  apparatus 
arrangement  is  shown  in  figure  1. 

Initial  testing  indicated  the  need  for  a  bench  alignment  procedure  for 
adjusting  the  periscope,  and  a  He-Ne  laser  alignment  system  was  developed. 
This  permits  overall  periscope  alignment,  including  optimizing  the  external 
adjustments  possible  at  the  last  external  Nd:YAG  beam  steering  mirror. 

Minor  adjustments  are*possible  with  the  periscope  system  installed  in  the 
vacuum  system  to  account  for  the  significant  differences  between  the  He-Ne 
and  Nd:YAG  laser  spot  sizes  and  beam  collimation. 


The  output  coupling  telescope  of  the  Nd:YAG  laser  has  been  adjusted  so 
that  the  frequency  doubled  Nd:YAG  beam  spot  size  is  about  one  cm  at  the  30 
cm  focal  length  lens  used  to  focus  the  laser  on  the  sample.  This  causes  the 
depth  of  focus  to  be  smaller  than  with  previous  experiments,  and  allows  a 
larger  laser  beam  spot  to  impact  the  periscope  optics,  with  less  likelihood 
of  laser  damage  to  the  optics  in  vacuum. 


The  final  element  of  the  periscope  is  a  high  laser  power  mirror, 
coated  for  532  nm  laser  radiation  at  45°  reflection.  Initial  experimental 
work  with  this  arrangement  indicated  that  vapor  deposition  and  laser  damage 
on  the  mirror  element  was  negligible.  However,  after  extended  use,  it  was 


noted  that  spot  damage  was  occuring,  but  with  Little  correlation  to 


conditions  of  operation  of  the  vacuum  svstem.  '.v’e  discovered  after  careful 
examination  of  the  laser  beam  distribution  that  the  damage  was  due  to  a 
laser  hot  spot.  It  appears  that  the  Nd:YAG  beam  had  been  developing  a  hot 
spot  for  some  period  of  time  which  coincidentially .  and  unfortunate Iv .  1-a 

to  a  best  tuning  pattern  similar  to  that  described  by  the  laser  operations 
manual.  The  damage  was  identified  as  a  series  of  small  burn  spots  on  the 
laser  rod.  This  resulted  in  an  inability  to  obtain  the  nearlv  gaussian  bea 
pattern  desirable  for  these  studies.  The  rod  and  a  damaged  polarizer  have 
been  replaced,  and  the  laser  system  is  again  fully  operational.  Also,  we 
plan  to  prevent  future  damage  by  using  a  coated  prism  as  the  final  element 


Modification  of  Signal  Detection  Apparatus 


As  part  of  the  new  experimental  arrangement,  we  have  implemented  a 
boxcar  integrator-based  scanning  technique  which  permits  a  slow  mass 
spectrometer  scan  of  an  entire  mass  spectral  region  synchronous  with  the 
laser.  Figure  2  shows  a  block  diagram  of  the  electronics  of  this 
experimental  arrangement.  The  ramp  generator  is  used  via  the  external 
command  input  of  the  quadrupole  control  to  sweep  the  mass  setting  of  the 
mass  spectrometer.  It  also  provides  the  trigger  to  start  the  sweep  of  the 
signal  averager  used  to  log  and  average  the  data  runs.  The  lamp  synch  from 
the  Nd:YAG  laser  triggers  the  time  base  of  the  boxcar  averager.  As  shown  i 
figure  2,  the  gate,  or  integration  time,  of  the  boxcar  is  positioned  in 
time  to  collect  the  ion  signal  synchronous  with  the  laser.  Since  the 
arrival  time  of  the  different  ions  are  mass  dependent,  the  width  of  the 
gate  has  to  be  set  wide  enough  to  incoporate  signals  for  all  masses  of 
interest.  This  unfortunately  increases  the  amount  of  integrated  background 


signal,  with  it's  inherent  noise,  relative  to  the  truly  synchronous  laser 
initiated  signal.  In  these  experiments,  with  the  preamp  input  resistance 
set  to  107  ohms,  the  gate  width  used  is  500  us.  At  20  Hz,  the  repetition 
rate  of  the  laser,  this  results  in  a  one  percent  duty  cycle  for  integrating 
the  signal.  The  dc .  non- laser  synchronous  component  could,  in  princple,  be 
removed  by  simultaneous  background  subtraction  using  the  second  gate  or 
channel  of  the  boxcar  averager.  However,  this  capability  was  not  used  since 
it  increases  the  inherent  ac  noise  in  the  signals. 

In  applying  a  deconvolution  model  to  the  time- resolved  data,  we 
previously  reported  that  the  mass  spectrometer  instrumental  response  was  a 
major  portion  of  the  convolution  function  and  would  need  to  be  dramatically 
improved  (at  least  by  two  orders  of  magnitude) ,  to  enable  extraction  of  the 
temperature-dependent  vapor-pressure  generation  function.  In  addition  to 
pressing  the  state  of  the  art  in  charge  amplifiers,  this  also  necessitated 
an  improvement  in  time  resolution.  We  have  addressed  this  problem  with  two 
approaches.  Initially,  the  mass  spectrometer  charge  amplifier  is  being 
operated  at  maximum  bandwidth  and  minimum  gain  (105  ohms).  Post¬ 
amplification  with  a  high  gain  differential  one  megahertz  bandwidth 
oscilloscope  amplifier  provides  the  high  level  signals  needed  by  the  signal 
and  boxcar  averagers.  We  are  currently  investigating  higher-bandwidth 
alternatives  to  the  existing  100%  feedback  charge  amplifier  technology.  The 
boxcar,  introduced  to  permit  gated  detection  while  sweeping  the  mass 
spectrum,  could  also  be  used  to  gather  high  speed  time - reso lved  data,  but 
only  at  a  very  low  duty  cycle,  i.e.  one  point  per  laser  shot.  To  improve 
time  resolution  we  have  initiated  a  procurement  for  a  much  faster  signal 
averager  to  match  the  higher  system  bandwidth.  As  an  alternative  to  analog 
data  acquisition,  we  have  also  contracted  for  a  custom  designed  ultra  high 
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speed  multichannel  scaling  system,  which  will  permit  ion  counting  with  time 
resolution  of  a  microsecond  or  less. 


In  the  new  perpendicular  sampling  geometry,  the  molecular  beam  pattern 
clearly  coincided  with  the  mass  spectrometer  beam  axis,  even  though  the 
periscoped  beam  was  directed  at  the  surface  at  approximately  a  15°  angle. 
Unlike  the  studies  at  15° ,  the  0°  studies  showed  significant  laser- 
synchronous  mass  spectrometer  peaks  at  the  background  mass  positions,  and 
at  mass  peaks  not  identified  with  carbon.  These  data  indicated  that  the 
carbon  Cn  (n=l-5)  peaks  were  of  similar  intensity  to  that  observed  off- 
axis,  even  though  estimates  of  macroscopic  carbon  transport  from  deposit 
patterns  indicated  more  than  an  order  of  magnitude  more  carbon  beam  densitv 
on  axis  relative  to  15°  off-axis.  Considerably  higher  than  normal  surface 
damage  was  also  noted,  which  was  subsequently  attributed  to  laser  hot- 
spotting  (from  the  defective  laser  rod). 


The  new  plume -samp ling  geometry  and  high  speed  signal  amplification 
has  allowed  us  to  monitor  the  evolution  of  surface  layer  and  impurity 
species  from  the  graphite  surface,  using  boxcar  integration  over  regions  of 
mass  spectra,  as  opposed  to  individual  mass  positions.  Three  typical  time¬ 
gated  mass  spectra  are  shown  in  figure  3.  In  figure  3a  the  laser  is  off  and 
the  mass  spectrometer  ionizer  is  on.  Thus  this  spectrum  shows  the 
background  signal  integrated  during  the  boxcar  gate.  The  relatively  high 
background  observed  here  is  abnormal,  and  was  identified  as  a  cryo  pumping 
problem  which  manifested  itself  by  an  inability  to  adsorb  permanent  gases 
on  the  time  scale  of  the  laser.  This  problem  has  been  resolved,  and  the 
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background  is  now  negligabLe  with  respect  to  the  plume  species.  The  second 
spectrum,  figure  3b,  has  the  laser  on  as  well  as  the  ionizer.  There  is  a 
laser- induced  increase  for  the  background  signals  noted  in  figure  3a.  In 
addition,  m/e  values  of  12,  23,  36  and  several  broad  peaks  above  m/e  of  50, 
appear  in  the  spectrum.  The  difference  between  these  spectra,  i.e.  figure 
3b  -  figure  3a,  would  be  the  laser  induced  mass  spectrum.  The  12  and  36  m/e 
values  are  identified  as  Cx  and  C3 ,  respectively.  These  spectra  serve  to 
illustrate  the  sizes  of  the  graphite  signals  relative  to  the  background.  In 
figure  3c  the  laser  is  left  on  but  the  ionizer  is  turned  off.  In  this  mode 
of  operation  the  signal  consists  only  of  ions  formed  synchronously  with  the 
laser.  Comparison  of  the  three  figures  indicates  that  the  signals  at  values 
m/e  of  23,  39,  and  41  are  mainly  the  result  of  ions  formed  during  the 
laser- induced  vaporization  process.  These  peaks  correspond  to  Na  and  the 
two  K  isotopes.  More  details  about  the  ion  detection  and  sensitivity  are 
discussed  below. 

In  the  case  of  K  at  m/e  39,  the  m/e  41  isotope  was  observed  with 
identical  peak  characteristics  to  m/e  39,  and  a  ratio,  I41/I39  of  0.08, 
essentially  identical  to  the  natural  abundance  ratio  of  0.074.  The  alkalies 
also  exhibited  the  bimodal  peak  shape  characteristic  of  other  non-carbon 
signals,  with  the  fast  component  being  dominant.  With  the  mass  spectrometer 
electron  bombardment  ion  source  off,  the  fast  component  of  Na  and  K  was 
identified  as  transmitted  ions  which  comprised  the  majority  of  the  mass 
spectrometer  signal.  A  small  portion  of  the  fast  component  of  the  carbon 
species  is  ions.  It  should  be  noted  that  the  mass  spectrometer  system  is 
intrinsically  about  four  orders  of  magnitude  more  sensitive  to  ions  than 
neutrals,  as  the  mass  spectrometer  electron  impact  ionization  process  is 
only  -0.01%  efficient.  These  ions  apparently  are  part  of  the  plume 


formation  process  and  should  serve  as  probes  of  the  plume.  Systematic 
measurement  of  ions  is  planned  when  we  have  the  special  ion  optics 
components  and  power  supplies  currently  in  procurement.  In  addition  to  the 
mass  spectrometer  studies,  time  -  resolved  optical  emission  (thermal)  and 
olume  fluorescence  studies  have  been  initiated. 

The  graphitic  carbon  signals  of  figure  4  show  distinctly  bimodal  peaks 
in  the  ion  intensity  vs.  time-of-arrival  data,  with  an  extremely  sharp 
initial  arrival  peak  which  we  were  unable  to  observe  with  the  slower  signal 
detection  electronics  used  previously  and  a  superimposed  slower  component. 
The  fast  component  of  these  signals  appears  to  be  the  result  of  ions  formed 
during  the  laser  vaporization  process  and  the  slower  component  is 
attributed  to  neutrals  formed  during  the  vaporization.  Ubiquitous 
background  peaks  now  also  appear  with  the  very  fast  component  and  a  broader 
component.  Additional  work  will  be  required  to  determine  how  much  of  the 
very  broad  component  of  these  surface  species  arises  from  electron  impact 
fragmentation,  or  other  effects. 

Results  of  laser  vaporization  on  graphite  samples,  which  were  polished 
to  ~20-25  micron  average  surface  roughness,  were  similar  to  those  for  the 
100  micron  surface  roughness  samples  studied  previously.  With  samples 
polished  to  the  order  of  a  few  microns,  we  can  begin  to  see  individual  shot 
damage,  as  expected.  We  experienced  some  sensitivity  to  beam  focus  in 
initiating  vaporization  with  low  angle  (i.e.  far  from  normal)  laser 
impingement  with  the  more  highly  polished  samples,  but  are  not  sure  how 
much  of  this  difficulty  resulted  from  the  erratic  beam  profile  behavior  due 
to  the  now-replaced  damaged  laser  rod.  Scanning  electron  microscope  work  to 
examine  the  microstructure  of  the  impact  zones  is  in  progress. 
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FIGURE  CAPTIONS 


Figure  1.  New  Experimental  Arrangement,  showing  beam  periscope  and  in¬ 
vacuum  x-v  translator  (dimensions  not  to  scale) . 

Figure  2.  Block  Diagram  of  Electronics  for  Scanned  Gated  Detection  of  Mas 
Spectra . 

Figure  3.  Scanned  Gated  Mass  Spectra  of  Graphite. 

Figure  4.  Time -Resolved  Mass  Spectrometer  Signals. 
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GATED  MASS  SPECTROMETER  SIGNALS 
(GATE  WIDTH  500  MICROSECONDS) 


